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I.  INTRODUCTION 


This  document  is  a  semiannual  report  prepared  under  Contract 
N00014-70-C-0161,  Modification  P00005.  The  purpose  of  the  program  per¬ 
formed  under  this  contract  is  to  continue  the  research  on  10.6  pm  image  up- 
conversion  detectors  and  to  determine  the  feasibility  of  developing  a  10.  6  pm 
range  gated  upconversion  imaging  laser  radar  receiver  with  150  by  150  reso¬ 
lution  elements.  To  meet  the  program  objectives,  the  following  study  areas 
have  been  required  by  the  contract: 

•  Experimental  determination  of  the  effect  of  a  diverging 
pump  on  image  resolution  for  the  image  space  configura¬ 
tion  upconverter  system,  and  comparison  of  the  results 
with  analytical  prediction. 

•  Analytical  determination  of  the  effect  of  birefringence 
and  a  truncated  pump  beam  on  image  resolution. 

•  Development  of  techniques  to  handle  power  dissipation 
inside  the  nonlinear  crystal. 

•  Investigation  of  techniques  for  improvement  of  optical 
efficiency. 

During  this  reporting  period  an  experimental  model  image  up¬ 
converter  capable  of  real  time  direct  viewing  has  been  assembled,  Imaging 
characteristics  such  as  resolution,  field  of  view,  and  the  upconversion  ef¬ 
ficiency  oi  the  experimental  model  were  measured  by  arranging  the  model  in 
the  Fourier  space  and  the  image  space  configurations. 

A  resolving  power  of  better  than  150  by  150  elements  has  been 
demonstrated  with  the  Fourier  space  configuration.  The  number  of  resolution 
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elements  obtained  with  the  image  space  configuration  is  65  by  65.  The  mea¬ 
sured  angular  resolution  in  both  configurations  is  nearly  equal  to  the  diffrac¬ 
tion  limit.  The  number  of  resolution  elements  obtained  with  the  image  space 
configuration  is  limited  by  the  area  of  usable  interaction  aperture. 

The  most  important  experimental  findings  with  regard  to  the 
imaging  properties  of  the  image  space  upconversion  are: 

•  Resolution  in  the  Fourier  space  configuration  is  degraded 
by  the  pump  beam  amplitude  distribution. 

•  Resolvable  image  linewidth  of  the  image  space  configura¬ 
tion  agrees  closely  with  the  prediction  based  on  the  phy¬ 
sical  optics  analysis.  The  geometric  optics  estimate  is 
three  to  four  times  poorer  than  the  measured  results. 

•  Resolution  in  the  image  space  configuration  would  not  be 
affected  by  the  pump  beam  intensity  profile.  An  effect 
of  the  pump  beam  intensity  distribution  appears  as  non- 
uniform  intensity  across  the  image  plane. 

•  Use  of  a  diverging  multimode  pump  beam  does  not  degrade 
image  resolution  if  the  upconverter  is  in  the  image  space 
configuration. 

•  Conversion  efficiency  of  the  TEMqo  mode  beam  pumped 

upconverter  approaches  the  estimate  (within  a  factor  of 
three)  for  perfectly  matched  plane  wave  interaction. 

The  experimental  model  was  also  operated  in  the  pulse  gated  mode.  Pre¬ 
liminary  results  of  this  experiment  are  also  presented  in  this  report. 

Analysis  describing  the  optical  transfer  function  of  the  Fourier 
space  configuration  and  the  line  spread  functions  of  non- Fourier  space  con¬ 
figuration  are  presented  in  Section  HI.  The  crystal  birefringency  is  included 
in  the  analysis.  Calculation  of  resolving  power  and  some  consideration  on 
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the  object  optics  are  also  included  in  Section  III.  The  main  results  of  this 
analysis  are: 

•  The  Fourier  space  configuration  image  upconverter  is 
indeed  free  of  aberrations.  The  field  of  view  is  limited 
by  the  phase  matching  condition  rather  than  the  crystal 
aperture  or  the  pump  beam  size. 

•  Resolving  power  of  the  Fourier  space  configuration  up¬ 
converter  pumped  by  a  Gaussian  beam  is  limited  by  the 
diffraction  due  to  the  Gaussian  aperture  when  the  lens 
diffraction  is  negligible. 

•  Because  of  the  crystal  birefringence,  the  upconverted 
image  of  a  line  object  normal  to  the  plane  containing  the 
crystal  C  axis  and  the  pump  wave  vector  consists  of 
several  narrow  nearly-equal-intensity  lines.  This  fine 
structure  has,  however,  not  been  observed  experimen¬ 
tally. 

Section  IV  presents  measurements  of  conversion  efficiency  and  realistic 
estimates  of  practically  achievable  efficiency.  It  is  concluded  that  the  con¬ 
version  efficiency  estimation  based  on  the  generally  accepted  plane  wave  in¬ 
teraction  could  be  achievable  at  the  center  of  the  field  of  view.  Realistic 
values  for  the  SHG  d  coefficients  of  proustite  are  also  considered. 

The  plans  and  the  areas  of  study  for  the  remainder  of  this  pro¬ 
gram  are  discussed  in  Section  V. 


II.  REAL  TIME  IMAGE  UPCONVERSION  EXPERIMENTS  WITH  THE 
DIRECT  VIEWING  LABORATORY  MODEL 

An  experimental  model  upconverter  capable  of  real  time  direct 
viewing  of  upconverted  images  has  been  set  up  during  this  period. 

The  imaging  characteristics  of  the  upconverter  were  measured 
first  with  the  proustite  crystal  pumped  by  a  CW  TEMqo  mode  beam.  The 
resolving  power  of  this  model  arranged  in  the  Fourier  space  configuration 
was  measured  to  assess  how  close  it  comes  to  achieving  diffraction  limited 
operation.  The  model  was  then  set  up  in  the  image  space  configuration  to 
measure  the  extent  of  thickness  aberration. 

The  measured  angular  resolution  of  approximately  1  milliradian 
for  the  Fourier  space  configuration  was  within  a  factor  of  two  of  the  Rayleigh 
resolution  of  a  uniformly  illuminated  object  lens  aperture.  The  number  of 
resolution  elements  contained  in  a  10-degree  field  of  view,  which  is  also 
experimentally  checked,  is  consequently  better  than  the  resolution  objective 
for  the  program  of  150  by  150  elements.  Upconverted  images  of  a  10.6-pm 
illuminated  USAF  1951  resolution  chart  containing  groups  smaller  than  1  lp/mm 
were  obtained.  Because  the  resolution  chart  does  not  fill  the  field  of  view, 
the  actual  number  of  resolvable  elements  contained  in  the  recorded  images 
which  fills  a  7. 5-degree  angle  is  114  by  114. 


The  resolving  power  of  the  model  arranged  in  the  image  space 
configuration  was  also  measured.  The  experimental  results  prove  (as 
expected)  that  the  resolvable  image  linewidth  is  limited  by  the  thickness 
aberration.  The  measured  image  resolution  linewidth  of  60  |im  agrees 
with  the  physical  optics  prediction  of  68  m,  but  is  approximately  three  to 
four  times  narrower  than  the  geometri'*  optics  calculation.  The  main  reason 
for  failure  of  the  geometric  optics  theory  to  accurately  predict  the  extent  of 
aberration  is  probably  due  to  the  inability  of  the  geometric  optics  to  account 
for  the  intensity  variation  of  the  upconverted  image  rays.  Upconverted  images 
of  the  resolution  chart  with  62  by  62  resolution  elements  were  obtained  with 
the  experimental  model  in  the  image  space  configuration.  Measured  angular 
resolution  is  slightly  worse  than  the  diffraction  liiwed  angular  resolution  of 
the  object  lens. 

The  effect  of  diverging  pump  beam  on  the  resolving  power  of  the 
upconverter  system  in  the  image  space  configuration  was  investigated.  The 
fact  that  resolution  of  the  image  space  configuration  system  using  a  multi- 
mode  pump  beam  does  not  degrade  from  that  of  a  TEMqq  mode  pumped  sys¬ 
tem  was  verified  experimentally.  However,  a  nonuniform  intensity  corre¬ 
sponding  to  the  multimode  pattern  is  superimposed  on  the  upconverted  image. 

Image  upconversion  experiments  with  the  image  space  configured 
model,  using  pulsed  pump  beam  and  a  gated  image  tube,  were  also  carried 
out.  Some  preliminary  results  of  the  experiments  are  presented  in  para¬ 
graph  E  of  this  section. 
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A.  EXPERIMENTAL  ARRANGEMENT 

The  experimental  model  upconverter  system  is  shown  schematically 
in  Figure  1.  The  system  was  arranged  so  that  the  optical  configuration  and 
the  pump  beam  parameters  could  be  readily  changed  with  minimum  physical 
alteration  of  the  setup.  A  2. 5-cm  diameter,  and  0. 6  cm  thick  proustite  crys¬ 
tal  was  used.  The  upconverter  was  aligned  for  Type  n  noncritical  phase 
matching.  The  proustite  was  oriented  such  that  its  crystallographic  X  axis 
lies  in  the  plane  of  paper  and  is  parallel  to  the  front  surface.  Its  optic  axis 
lies  in  a  plane  perpendicular  to  the  paper  and  makes  an  angle  of  20  degrees 
with  the  normal  (along  the  C  axis)  to  the  polished  front  and  back  surfaces. 

A  continuously  pumped  Nd:YAG  laser  radiating  electric  fields,  polarized 
parallel  to  the  plane  of  paper,  was  used  as  a  pump  source.  The  laser  is 
fitted  with  an  acousto -optic  Q-switch  allowing  pulsed  as  well  as  CW  opera¬ 
tion.  The  beam  parameters  of  the  1. 06-jx  m  pump  beam  emanating  out  of  the 
team  expander  are  shown  in  Table  1  for  various  modes  of  operation. 

TABLE  1.  1.06 -um  PUMP  BEAM  PARAMETERS 
AFTER  BEAM  EXPANSION 


MODE 

AVERAGE 

POWER 

BEAM 

CROSS  SECTION 

BEAM 

DIAMETER 

BEAM 

DIVERGENCE 

CW  TEMqq 

5.0  watts 

Gaussian 

1. 1  cm 

0. 16  mrad 

CW  Multi¬ 
mode 

7. 5  watts 

1. 1  cm 

1.0  mrad 

Q-Switched 

TEMoo 

1000  PPS 

0. 5  watt 

Gaussian 

1. 1  cm 

0. 16  mrad 
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A  polarizer  external  to  the  NdtYAG  laser  removes  field  com¬ 
ponents  polarized  along  the  y  direction.  The  chopper  reduces  the  average 
power  of  the  pump  beam  by  serving  as  an  attenuator  when  required.  Glass 
lenses  and  Lg  collimate  and  expand  the  pump  beam. 

A  multimode  CC>2  laser  radiating  with  an  electric  field  polarized 
perpendicular  to  the  plane  of  the  paper  was  used  to  back  illuminate  the  trans¬ 
parent  1951  USAF  resolution  chart.  The  COg  beam  diameter  at  the  resolu¬ 
tion  chart  is  sufficiently  large  to  illuminate  the  entire  resolution  chart  by 
adjusting  the  10.6-pm  path  length  by  means  of  mirrors  as  shown  in  Figure  1. 
The  position  of  BaFg  lens  L^,  relative  to  the  proustite  crystal  and  the  object, 
is  determined  by  the  configuration  of  the  system  (image  space  or  Fourier 
space). 

The  pump  wave  and  the  object  wave  are  combined  by  a  dichroic 
mirror.  An  improved  dichroic  mirror  was  obtained  during  this  reporting 
period.  This  mirror  reacts  more  than  99  percent  of  the  YAG  beam  and  has 
a  transmission  loss  at  10. 6  pm  of  about  13  percent.  As  shown  in  Figure  1, 
the  angle  of  incidence  of  the  YAG  beam  on  the  dichroic  mirror  is  about  20°. 
This  eliminates  potential  vignetting  of  10. 6-pm  object  waves  by  the  mirror 
holder. 

A  Gian  air  polarization  filter  and  narrowband  interference  filters 
are  used  for  suppression  of  1. 06-pm  pump  wave.  Transmission  through 
these  filters  at  0.967  pm  is  about  10  percent.  Lens  is  used  to  form 
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upconverted  real  images  on  the  photo  surface  of  the  image  intensifier  tube. 
The  image  intensifier  tube  is  an  ITT  Model  4740  three-stage  gated  tube  with 
an  S-l  photocathode.  Its  pertinent  characteristics  are  listed  in  Table  2. 

The  image  intensifier  tube  can  be  cooled  to  -30°C  with  dry  ice  to  reduce 
tube  thermal  noise. 


TABLE  2.  CHARACTERISTICS  OF  ITT  MODEL  F4740  IMAGE  TUBE 


Useful  Cathode  Diameter 
Radiant  Sensitivity  at  0. 967  pm 
Center  Resolution 
Peripheral  Resolution 
Paraxial  Magnification 
Distortion 
Luminous  Gain 

Conversion  Efficiency  at  1.06  pm 

Radiant  Power  Gain  Peak  Wavelength 

Radiant  Power  Gain  1.06j* 

Equivalent  Background  Energy  Input  at 
Peak  Wavelength 

Equivalent  Background  Energy  Input  at 
1.06/x 


23  mm 
1.3  mA/W 
28  lp/mm 

25  lp/mm  at  7.0-mm  radius 
0.84 

18  percent  at  10-mm  radius 
6000  fL/fc 

1  000  ooo  — ^~2^)  ■ 

i,ui  u?uuu  W(106) 

1040  W/W 
208  W/W 

2  x  10“9  W/cm2 

8  x  IQ"9  W/cm2 


To  check  the  resolving  power  of  the  0.967  gm  image  detector 
used  in  the  system,  consisting  of  lens  L^,  the  interference  filter,  the  image 
intensifier  tube,  and  the  camera,  the  resolution  chart  is  illuminated  with 
visible  light  and  imaged  on  the  photosurface  of  the  image  intensifier  tube. 
Figure  2  shows  a  typical  image  recorded  on  the  film.  The  smallest 
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RESOLUTION  LIMITS  OF  THE  0.967-ym  IMAGE 
DETECTOR  OF  THE  EXPERIMENTAL  MODEL  IS  37  ym 
AS  SEEN  IN  THE  ILLUSTRATION 


FIGURE  2.  IMAGE  OF  THE  RESOLUTION  CHART  ILLUMINATED  BY 
VISIBLE  LIGHT 


resolvable  linewidth  on  the  output  phosphor  of  the  image  tube  is  37  yin. 

A  modified  polaroid  scope  camera  was  used  to  record  the  upconverted  images 
produced  by  the  image  tube. 


B.  IMAGE  UPCONVE RSION  IN  FOURIER  SPACE 

The  limiting  resolving  power  of  the  experimental  model  was 
measured  by  operating  the  model  in  Fourier  Space  configuration  and  pump¬ 
ing  the  nonlinear  material  with  a  TEM  mode  beam.  Diffraction  limited 
^  oo 

image  upconversion  is  possible  for  this  configuration  if  the  pump  is  a  plane 
wave^.  Thickness  aberrations  are  negligible  for  this  configuration. 

A  1951  USAF  resolution  chart  was  placed  in  the  front  focal  plane 
of  lens  Lg.  The  distance  between  BaF2  lens  Lg  and  the  proustite  crystal  is 
120  mm  which  is  equal  to  the  focal  length  of  lens  Lg.  Lens  is  placed  so 
that  the  proustite  crystal  and  the  photocathode  of  tne  image  intensifier  tube 
coincide  with  the  front  and  back  focal  pianes  of  lens  L^,  respectively.  The 
focal  length  of  lens  is  580  mm. 

The  pump  beam  diameter  at  the  proustite  front  face  is  1. 1  cm  and  its 
beam  divergence  is  0. 16  mrad.  The  pump  beam  is  in  TEMqq  transverse  mode. 
Pictures  reproduced  from  some  of  the  upconverted  images  and  recorded 
by  the  polaroid  camera  are  shown  in  Figure  3.  The  smallest  element  that 
can  be  visually  resolved  in  these  figures  is  the  group  2/1.  Figure  3A  has 
better  resolution  along  the  x  direction,  but  Figure  3B  has  better  resolution 
along  the  y  direction.  (See  Figure  *  for  the  coordinate  system.)  Inspection 


of  upconverted  images  on  the  output  phosphor  surface  of  the  image  intensifier 
tube  indicates  that  group  1/6  is  resolved  along  both  the  x  and  ihe  y  directions. 
The  system  angular  resolution  is  calculated  by  the  formula 


y  =  h.  /fQ 
r  it  3 


(2-1) 


where 


h.  =  resolvable  linewidth  of  the  test,  chart 
it 


fg  =  effective  focal  length  of  lens 


The  measured  angular  resolution  corresponding  to  the  1/6  and  2/1  resolu¬ 
tion  groups  is: 


*-3 


Group  1/6  h^  =  1/7. 13  mm  =  1. 17  x  10"°  rad 


Grouo  2/1  h.  -  1/8  tnm  =  1.04  x  10~3  rad 


IT 


Measured  angular  resolution  of  1. 17  mrad  is  2.  5  times  larger  than  the 
Rayleigh  resolution  angle  of  3-cm  diameter  lens  L^. 

Resolution  degradation  from  the  Rayleigh  resolution  angle  of  the 
object  lens  is  mainly  due  to  nonuniform  transverse  intensity  distribution  of 
the  pump  beam.  It  will  be  shown  later  that  a  Gaussian  profile  of  the  pump 
beam  would  introduce  an  additional  Gaussian  attenuation  factor  on  the  angular 
spectrum  of  the  upconverted  image  field.  It  is,  therefore,  expected  that  the 
resolving  power  of  the  experimental  model  is  limited  by  the  Gaussian  window 
diffraction  rather  than  the  object  lens  diffraction. 
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The  system  field  of  view  was  measured  by  translating  the  object 

along  the  x  and  y  directions  until  the  edge  of  the  test  chart  disappeared. 

The  measured  system  FOV  of  0. 183  rad  (10. 5°)  agrees  with  the  previous 
(1  2  3) 

result'  *  *  .  Figure  3A  shows  that  the  FOV  covers  more  than  the  area 

occupied  by  the  resolution  chart.  The  number  of  resolvable  elements  con¬ 
tained  in  the  FOV  is  found  to  be  156  by  156,  if  one  takes  group  1/6  as  the 
minimum  resolvable  element.  The  image  shown  in  Figure  3A  contains 
approximately  114  by  114  resolution  elements. 

The  measured  transverse  magnification  of  0.99  agrees  within 
10  percent  of  the  estimated  value  of  1.08,  which  is  calculated  by  the  ex¬ 
pression 

MHM«cMcMt| 

where 

*s  f4  \ 

M  =  - —  -=-=0.44  :  upconversivn  magnifies :iov. 
uc  Air  3  / 

M  (=2.94):  camera  magnification 
c 

M,j,  (  =  0. 81):  image  tube  magnification 

C.  IMAGE  UPCONVERSION  IN  IMAGE  SPACE 

Image  upconversion  experiments  were  carried  out  with  the 
experimental  model  arranged  in  the  image  space  configuration.  Lens  of 
Figure  1  was  positioned  such  that  a  real  10. 6-^m  image  of  object  is  focused 
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onto  the  crystal.  The  photosurface  of  the  image  intensifier  and  the  proustite 
crystal  coincide  with  the  conjugate  planes  of  lens  L^. 

It  has  been  shown  analytically  that  the  resolvable  image  linewidth 

of  the  image  space  configuration  upconverter  pumped  with  a  plan*  wave  is 

/i  0} 

limited  by  thickness  aberration^  ’  .  It  is  the  purpose  of  this  experiment  to 

measure  the  aberrations  and  to  determine  the  limiting  resolution  of  the  image 
space  configured  upconverter.  The  Nd:YAG  laser  was  operated  in  the  TEMqo 
mode  to  minimize  the  effect  of  diverging  pump  beam. 

Figure  4  shows  reproduced  pictures  of  upconverted  images  re¬ 
corded  by  the  camera.  The  magnification  of  image  forming  lenses  and 
are,  respectively,  0.304  and  2.76.  The  smallest  resolvable  group  (1/3)  has 
a  linewidth  of  2. 52  lp/mm.  The  corresponding  image  width,  calculated  by 
dividing  the  resolvable  target  width  by  the  magnification  factor  of  lens  L^,  is 
60.3  pm. 

An  experimentally  measured  resolvable  image  linewidth  of  60  pm 
is  approximately  a  factor  of  3  narrower  than  the  resolvable  image  width  pre¬ 
dicted  by  the  geometric  optics  calculation.  As  will  be  shown  in  Section  m, 
the  physical  optics  calculation  predicts  that  the  minimum  resolvable  image 
width  for  the  given  parameter  values  is  68  pm.  The  experimental  result 
agrees  better  with  the  physical  optics  prediction  than  with  the  geometric , 
optics  calculation. 


FIGURE  4.  UPCONVERTED  IMAGE  OF  1951  USAF  TEST  CHART 
PRODUCED  BY  IMAGE  SPACE  UPCONVERTER  WITH 
SINGLE  MODE  PUMP 


In  the  geometric  optics  model,  the  analysis  is  based  on  the 
assumption  that  each  infinitesimal  slab  of  material  is  capable  of  forming  an 
independent  upconverted  image.  The  concept  of  thickness  aberrations  is, 
therefore,  introduced  by  considering  the  upconverted  image  as  a  congregation 
of  images  formed  by  infinitesimal  slabs  of  upconverter  material.  It  has  been 
previously  shown  (by  the  Fourier  optics  treatment)  that  the  aberration  is 
physically  introduced  by  imperfect  phase  matching  of  angular  spectra  of  the 
interacting  electric  fields.  The  angular  spectrum  with  higher  spatial  fre¬ 
quency  suffers  higher  attenuation  due  to  the  phase  mismatch.  The  degree  of 
phase  mismatch  is  proportional  to  the  material  thickness.  In  the  physical 
optics  model  the  upconverted  image  is  formed  by  linear  superposition  of  the 
upconverted  rays  (angular  spectra)  emanating  from  the  exit  surface  of  the 
material.  It  is  expected  that  the  aberrations  calculated  by  these  two  models 
should  be  different  from  each  other.  The  experimental  results  should  agree 
better  with  the  Fourier  optics  predictions  rather  than  the  geometric  optics 
because  the  model  is  in  better  agreement  with  the  actual  physical  process. 

The  number  of  resolution  elements  contained  in  the  image  shown 
in  Figure  4  is  approximately  62  by  62  which  is  approximately  a  factor  of  two 
less  than  that  obtained  by  the  Fourier  space  configuration  system.  Measured 
system  magnification  is  approximately  10  percent  greater  than  the  theoretical 
estimate. 

The  system  angular  resolution  can  be  calculated  by  dividing  the 
measured  resolvable  image  width  by  the  distance  between  the  crystal  and 
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lens  Lg.  Measured  system  angular  resolution  of  0. 4  mrad  is  very  close  to 
the  Rayleigh  resolution  angle  of  0.42  mrad  of  lens  L^.  The  experimental 
model  is,  therefore,  nearly  diffraction  limited. 

It  must  be  noted  that  the  Gaussian  distribution  of  the  pump  beam 
does  not  degrade  the  resolution  as  expected  from  the  physical  optics  analysis. 

It  can  easily  be  shown  that  the  Gaussian  amplitude  factor  of  the  angular 
spectrum  of  the  upconverted  image  fields  becomes  unity  for  all  angular 
frequencies  in  the  image  space  configuration  where  the  apparent  IR  object 
location  coincides  with  the  nonlinear  material^ 

The  primary  effect  of  Gaussian  profile  of  the  pump  beam  appears 
as  intensity  reduction  of  upconverted  images  toward  the  edge  of  FOV.  The 
size  of  the  10. 6  -p  m  image  of  the  object  formed  at  the  crystal  would,  there¬ 
fore,  be  considerably  smaller  than  the  pump  beam  size  if  one  wants  to  avoid 
intensity  reduction  at  the  outer  edge  of  FOV.  In  the  experiment,  the  size  of 
the  10, 6 -pm  image  of  the  resolution  chart  at  the  crystal  was  3.6  mm  by  3. 6  mm 
(considers iny  smaller  than  the  pump  beam  diameter).  Loss  of  brightness  and 
resolution  at  the  edge  can  even  be  seen  in  Figure  4,  which  was  taken  with 
some  precaution  to  minimize  the  Gaussian  window  effect.  If  the  pump  beam 
is  expanded  to  1-inch  diameter  and  usable  pump  power  is  increased  by  four 
times  that  of  the  present  system,  the  number  of  resolution  elements  could 
be  increased  to  a  minimum  of  130  by  130,  since  the  resolvable  image  width 
remains  the  same. 


< 


! 
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It  was  proved  experimentally  that  the  <  *ble  image  width 
was  unchanged  while  the  magnification  of  lens  did  change. 

D.  IMAGE  SPACE  UPCONVERSION  WITH  A  MULTIMODE  PUMP 
BEAM _ 

The  purpose  of  this  experiment  is  to  prove  that  a  multimode 

pump  beam  can  be  used  in  the  image  space  configuration  system  without 

sacrificing  the  resolving  power  from  the  single  mode  pumped  image  space 

system.  It  is  widely  recognized  that  the  image  resolution  degrades  with  an 

(o  4) 

increase  of  pump  beam  divergence  ’  .  This  has  been  demonstrated 

experimentally  for  the  Fourier  space  configuration  upconverter^  and  for 
the  optical  system  in  which  an  apparent  IR  object  is  located  at  finite  distance 
from  the  nonlinear  crystal^’ 

The  reduction  of  image  resolution  due  to  the  pump  beam  diver¬ 
gence  (or  equivalently  the  pump  beam  profile)  has  been  analyzed  by  the  ray 
tracing  technique^1  ’  ^  and  the  physical  optics  approach^.  It  has  been 
shown  that  the  amplitude  of  the  angular  spectrum  of  the  upconverted  image 
of  an  on  axis  point  object  contains  an  attenuation  factor  that  depends  on  the 
pump  beam  intensity  distribution  and  the  object  position^.  This  attenuation 
factor  becomes  independent  of  the  pump  beam  profile  when  the  apparent  IR 
object  plane  coincides  with  the  equivalent  pump  beam  plane.  The  resolving 
power  then  becomes  insensitive  to  the  pump  beam  profile  and  will  be  limited 
by  the  thickness  aberration.  It  must  be  understood,  however,  that  the 
intensity  distribution  of  the  upconverted  image  depends  strongly  on  the  pump 
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beam  intensity  distribution.  In  order  to  prevent  undesirable  image  intensity 
variation  across  the  FOV  in  the  image  space  configuration  system,  a  ran¬ 
domly  phased  multimode  pump  beam  is  preferred  over  a  beam  in  a  single 
(or  a  few)  higher  order  modes  which  may  have  holes  in  the  beam. 

Figure  5  shows  reproduced  pictures  of  an  upconverted  image  of 
the  resolution  chart.  The  image  shown  in  Figure  5A  was  obtained  with  a 
TEMqo  mode  pumped  image  space  system.  Figure  5B  shows  an  image  ob¬ 
tained  with  the  same  system  pumped  with  a  multimode  beam.  The  two  images 
remain  essentially  the  same,  indicating  no  degradation  of  resolution  due  to 
the  pump  beam  divergence.  It  was  observed,  however,  that  the  image  ob¬ 
tained  from  the  multimode  pump  beam  was  less  bright  than  that  produced  by 
the  TEMqo  pumped  beam.  Quantitative  comparison  of  image  brightness  for 
these  two  cases  was  not  carried  out. 

It  appears  that  dependence  of  upconversion  efficiency  upon  the 
pump  beam  divergence  requires  further  analytical  and  experimental  investiga¬ 
tion  to  design  an  efficient  image  upconversion  system. 

E.  IMAGE  UPCONVERSION  WITH  A  PULSED  PUMP  LASER  AND  A 
SYNCHRONOUSLY  GATED  IMAGE  TUBE _ 

An  experimental  model  pulse-gated  image  upconversion  system 
was  assembled,  and  some  preliminary  experiments  were  carried  out.  The 
optical  layout  of  the  model  is  identical  to  the  image  space  configuration  (Fig¬ 
ure  1).  A  block  diagram  of  electronics  circuits  used  for  the  synchronous 
pump  pulse  generation  and  gating  of  the  image  tube  is  shown  in  Figure  6. 


Nd:  YAG  LASER 


TO  IMAGE  INTENSIFIER 
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FIGURE  6.  BLOCK  DIAGRAM  OF  SYNCHRONOUS  GATE  ELECTRONICS 


The  Nd:YAG  laser  operating  in  single  transverse  mode  was  Q-switched  by 
an  acousto -optical  Q  switch.  The  pulse  repetition  rate  is  controlled  by  the 
trigger  pulse  generator  and  can  be  varied  from  1  to  5  kHz.  Pump  pulses 
with  200-ns  duration  and  approximately  5-kW  peak  power  were  expanded  to 
1-cm  diameter.  No  damage  to  the  crystal  was  observed.  The  high  voltage 
gate  pulse  generator  has  a  40-ns  rise  time  and  an  80-ns  fall  time.  Puise  width 
and  the  pulse  position  of  the  gate  pulses  can  be  adjusted  for  the  optimum 
image  signal-to -noise  ratio. 

The  images  obtained  with  the  synchronously  gated  image  tube 
are  shown  in  Figure  7A.  The  image  shown  in  Figure  7B  is  obtained  by  the 
same  system,  but  the  image  tube  was  not  gated.  The  improvement  in  the 
image  signal-to -noise  ratio  attained  by  gating  the  image  is  quite  remarkable. 
Optimization  of  the  experimental  system  and  quantitative  measurement  of 
image  resolution,  conversion  efficiency,  image  SNR,  etc.,  has  begun. 


IH.  ANALYTICAL  CONSIDERATIONS 


The  imaging  properties  of  the  parametric  image  upconversion 
process  have  been  previously  analyzed  by  the  geometric  optics  (ray  tracing) 
technique  *  ,  and  the  Fourier  optics  technique  v  '.  It  was  shown  that  the 
resolving  power  of  the  image  upconverter  is  (neglecting  the  aperture  diffrac¬ 
tion)  limited  by  thickness  aberration  and  thickness  coma.  The  aberrations 
are  caused  by  the  phase  match  error  inherent  in  the  wide  angular  spectral 
of  object  waves.  Analytical  expressions  for  the  calculation  of  abberations, 
and  equations  relating  the  angular  spectral  of  the  upconverted  image  field 

with  object  wave  Fourier  components  for  different  upconverter  configura- 

(2) 

tions,  were  previously  reported 

It  is  the  purpose  of  this  section  to  extend  previous  analysis  for 
further  understanding  of  the  parametric  image  upconversion  process.  The 
point  spread  function  and  the  optical  transfer  function  of  a  Fourier  space 
configuration  upconverter  system  is  discussed  first.  The  result  confirms 
an  earlier  conclusion  of  aberration  free  operation  with  the  Fourier  space 
configuration.  Considerations  of  the  resolving  power  of  the  image  space 
configuration  upconverter  systems  with  plane  wave  pump  are  presented  next. 
This  treatment  follows  the  geometric  optics  approach.  The  effect  of  mat¬ 
erial  birefringence  is  considered  last.  Line  spread  functions  of  the  plane 
wave  pumped  upconverter  system  are  obtained  without  neglecting  the  bire¬ 
fringence.  Image  blurring  due  to  the  crystal  birefringence  could  be  non- 
negligible. 
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A.  POINT  SPREAD  FUNCTIONS  AND  OPTICAL  TRANSFER  FUNCTION 
OF  A  PLANE  WAVE  PUMPED  INFINITE  APERTURE  FOURIER  CON¬ 
FIGURATION  SYSTEM 


Diffraction  due  to  apertures  is  neglected  in  this  discussion  in 
order  to  separate  image  properties  of  the  upconversion  process  from  the 
diffraction.  Lens  apertures  (the  nonlinear  material  aperture)  and  the  pump 
beam  diameter  are,  therefore,  assumed  to  be  infinitely  large. 

A  simplified  schematic  of  a  Fourier  space  configuration  upcon- 
verter  system  is  shown  in  Figure  8.  Lens  L^  is  positioned  such  that  the 
object  and  the  front  surface  of  the  nonlinear  material  are  respectively  in 
the  front  and  back  focal  planes.  The  front  local  plane  of  lens  L^  coincide 
with  the  exit  surface  of  the  nonlinear  crystal.  Upconverted  real  image  then 
forms  in  the  back  focal  plane  of  lens  Lg.  Lenses  Lj  and  L2  therefore, 
form  a  Fourier  transform  pair.  The  upconversion  takes  place  in  the 
Fourier  transformed  space  of  the  object  waves;  thus  the  name  Fourier  space 
configuration  system. 

For  the  plane  wave  pumped  upconverter,  it  can  be  shown  that 
the  angular  spectrum  of  upconverted  image  fields  in  air  beyond  the  non¬ 
linear  crystal  (z  >  +  D/2)  is  related  to  the  Fourier  component  of  object 

(1) 


fields  over  the  front  surface  of  the  crystal  by 


n 


Ug,  z)  ~  £.  (*s,-D/2)  Sine  ^ 
elPszD  +  Ksw  <z  "  °/2D 


-i  AKD/2 


(1) 
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FIGURE  8.  IMAGE  UPCONVERTER  ARRANGED  IN  FOURIER  SPACE 
CONFIGURATION 
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whore 


Sine  x  =  (yin  x)/x 
D  :  material  thickness 


z  component  of  wave  vector  at  X  in  the 
material  s 


K  _  :  z  component  of  wave  vector  at  X  in  air 

SoZ  S 


A  K  (=  Ksz  -  K. ,  -  t  k  v  +  _ )  :  phase  mismatch  factor 

2K-, 

M 

Km  :  wave  vector  at  Xg  along  the  phase  matched  direction 
C  :  anisotropy  angle 

The  angular  spectrum  £.  (  k  ,  -  D/2)  is  a  Fourier  transform  of  the  object 
field  distribution  over  the  plane  z  =  -D/2.  The  object  field  distribution  in 
the  plane  z  =  -D/2  is  a  Fourier  transform  of  the  object  field  distribution  over 
the  front  focal  plane  of  lens  since  the  plane  z=-D/2  coincides  with  the 


back  focal  plane  of  lens  Lj  '  . 


r  -i  j  p  .  _p.  j 

Ei  (P ’  'D/2)  =  r"  Ei  (gp  "2fl  '  D/2>  e  1  dii  (2) 


fii  (ft  -  °/2)  =  7j 


E.  (p,  -  D/2)  e 


-1  K  •  p 


dp  (3) 


. . . 


where 


f  j  =  focal  length  of  lens 

jp,  =  transverse  position  in  the  front  focal 
plane  of  lens 

p  =  transverse  position  in  the  back  focal  plane 
(z  = -D/2)  of  lens 

Substituting  Equation  (2)  into  Equation  (3)  and  carrying  out  the 
integration  we  obtain 

£.(*:,  -D/2)~E.  (-f2  /K.a,  -2f t  -  D/2)  (4) 

Constant  Amplitude  factor  and  fixed  phase  factors  have  been 
omitted  in  Equation  (4)  since  they  do  not  affect  the  result  in  any  significant 
way.  Equation  (4)  indicates  that  the  angular  spectrum  of  field  distribution 
over  the  back  focal  plane  of  a  nondiffracting  thin  lens  has  the  same  functional 
form  as  the  field  distribution  over  the  front  focal  plane  of  the  lens.  This 
result  is  no  surprise  if  one  realizes  that  they  are  related  by  double  Fourier 
transform. 

The  field  distribution  of  an  upconverted  image  formed  in  the 
back  focal  plane  of  lens  L0  has  the  same  functional  form  as  £  (  £  ,  -D/2) 
since  the  plane  z  =  +D/2  and  the  upconverted  image  plane  are  two  focal 


planes  of  lens  Lg. 


(5) 


Esks’2f2  +  D/2>~  4(Ksa<>s/f2’D/2> 

Combining  Equations  (1),  (4),  and  (5),  we  obtain  an  expression  for  the 
upconverted  image  field  distribution 


K  f, 
sa  1 


E  (p  2  f  +  D/2)  ~  E  -  2L  -  D/2). 

S  _s  1  1  Kiaf2  1  (6) 


Sine  (AP)  e 


-l 


AKD 


l  2 


-  K  D 
sz 


where  mismatch  factor  AK  and  the  z  component  wave  vector  K  in 

SZ 


Equation  (6)  are  evaluated  by  replacing  *  with  K  „  p  /f0  and  are  given  by 

S  Sa — S  u 

2 

i  i 

(V 


K 


AK  =  Ks  -  -  C  Ksa  xs/f2  + 


Sa  -  „  2  (2  .1) 
2  f  2  s  KM  Ks 


and 


Ksz2-Ks2-<Ks^/f2>2 


(8) 


If  crystal  anisotropy  is  neglected  the  mismatch  factor  aK  is 


K 


AK  =  - 


The  function  sine 
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sa 


\f2 


(9) 


’  K  D 

P 

,Ksa  1 

2 

[4KMKi 

1*2  P-J 

- 

indicates  that  the  amplitude 


of  an  upconverted  image  field  of  a  uniformly  illuminated  object  field  will  not 

be  uniform  but  will  decrease  with  transverse  distance  p  .  It  imposes  a 

s 

limitation  on  the  size  of  the  upconverted  image  but  will  not  limit  the  image 
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resolution.  If  we  define  the  size  of  an  upconverted  image  as  equal  to  the 

transverse  distance  p  ,  where  the  first  zero  of  a  sine  x  function  occurs, 

s 

the  radius  of  an  upconverted  image  is  expressed  by 


K. 

1 


(10) 


Equation  (6)  also  shows  that  the  transverse  magnification  of  a  Fourier  space 
configuration  upconverter  system  is 


Kiaf2 

KsafX 


Xs  f2 

vr 


(id 


Equation  (11)  can  also  be  obtained  by  the  geometric  optics  analysis.  The 
limiting  size  of  an  object  that  can  be  upconverted  can  be  easily  found  by 
combining  Equations  (10)  and  (11). 

f,  InK^K. 

VVN  =2  <12> 

ia  y  P 


Noting  that  the  object  size  a.  is  related  to  the  acceptance  angle  0.  (FOV)  of 

1  13. 

the  upconverter  by 


we  obtain 


/  2  A  X.n.n 
p  l  i  sM 

1  AnD 
\  s  P 


1/2 


(13) 
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For  a  noncritically  phased  matched  proustite  unconversion  from  10. 6  pun  to 

1  11/2  H 

D  rad' 


0. 967  urn,  we  have  0.  =  0. 157 

la 


To  illustrate  the  points  presented  in  the  previous  discussion  let  us  consider 
a  case  where  the  object  field  distribution  is  a  sine  function  as  shown  in 
Figure  9A.  The  field  distribution  of  upconverted  image  formed  in  the  back 
focal  plane  of  lens  L2  is  shown  in  Figure  9B.  The  period  of  the  sine  function 
image  field  is  magnified  by  a  factor  and  its  amplitude  envelope  follows 
the  sine  function.  The  image  is  inverted  as  shown  in  Figure  9B  so  that  an 
object  point  -y^  is  imaged  at  ygl  =  y.^.  The  resolving  power  is  un¬ 

limited  since  the  image  field  period  is  linearly  proportional  to  the  object 
field  period.  The  object  size  that  can  be  upconverted  is,  however,  limited. 
The  size  limitation  is  caused  by  the  phase  match  requirement  rather  than 
the  crystal  aperture  limitation. 

The  point  spread  function,  which  is  a  response  to  a  point  object, 
is  easily  obtained  from  Equation  (6).  The  field  distribution  of  a  point  object 
at  (x.q,  y.o)  is  expressed 

E,  (Pt)  ~  6  (Xj  -  xi0)  •  6(Yi  -  yto) 


Then  from  Equation  (6),  the  point  spread  function  becomes 

Hs  fes*  ~ 6  <VMt  -  xio>  6(ys/Mt  -  yJ 


-i  (AK/2  -  Kgz)  D 


Sine 


(14) 


J*=KsaMtAio/f2 
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A.  OBJECT  FIELD  DISTRIBUTION 


B.  UPCONVERTED  IMAGE  FIELD  DISTRIBUTION 

2-4815 

FIGURE  9.  FIELD  DISTRIBUTION  OF  OBJECT  AND  UPCONVERTED 
IMAGE  FOR  THE  FOURIER  SPACE  CONFIGURATION 


where 


—  =  Ksa£s  ^2  in<iicates  <k^at  the  Quantity  inside  the 
bracket  is  evaluated  at  *.  =  K  p  /f  as  shown  in  Equation  (7)  and  Equation  (8). 

Set  —  S  u 

The  coherent  transfer  function,  which  is  the  Fourier  transform 
of  the  point  spread  function,  becomes 


„  ,  .  m2  -1Mt^s  'fio 

hs(is)~Mt  e 


Sine 


AKD 
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(4£-k  ] 
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2  sz 

(15) 


k  =K  M ,p.  /U 

sa  t  —io  2 


Amplitude  of  h  (  k  )  depends  on  the  object  location  and  thus  object  size 
limitation  is  implied  even  for  the  infinite  aperture  Fourier  configuration 
upconverter.  The  optical  transfer  function,  which  is  defined  ^  as, 

fs  (*.')•  h*  U' +  *  ><*J. 


W3— 5 


//»s  («' 


)  *  d  k' 


becomes 


Hs  <£.)  =  e 


i  M,  k  •  p. 
t  _  —10 


(16) 


(17) 


Equation  (17)  again  indicates  the  property  of  unlimited  resolving  power  since 
the  amplitude  of  optical  transfer  function  (Modulation  transfer  function)  is 
independent  of  angular  frequency  k  . 


B.  RESOLVING  POWER  OF  THE  IMAGE  UPCONVERTER  IN  NON-FOURIER 
SPACE  CONFIGURATION _ _ 

Resolving  power  of  the  plane  wave  pumped  image  upconverter 

systems  is  considered  in  these  paragraphs.  The  crystal  birefringence  is 

neglected  and  paraxial  assumption  is  used.  Computational  procedure  based 

on  the  ray  tracing  approach  has  been  suggested  previously  without  detailed 
.  (2) 

discussions  .  The  present  consideration  is  intended  to  give  computational 
formulas  for  various  optical  arrangements. 

(2) 

It  has  been  shown,  both  from  the  geometric  optics  and  the 
physical  optics  considerations  ^ ,  that  the  minimum  resolvable  image  width 
of  the  non-Fourier  space  configuration  is  independent  of  the  longitudinal 
object  location  if  the  thickness  aberration  is  the  only  factor  affecting  image 
resolution.  It  has  also  been  shown  that  transverse  magnification  of  the 
plane  wave  pumped  upconverter  is  independent  of  object  location  and  is  equal 
to  the  resolvable  upconverted  image  width. 

Resolution  of  two  types  of  non-Fourier  space  configuration  is 
considered  based  on  the  geometric  optics  approach.  Two  optical  systems, 
following  Andrew’s  designation,  are  schematically  shown  in  Figure  10.  The 
optical  system  n  (Figure  10A)  has  the  following  arrangement. 

A  real  object  is  located  a  distance  S  from  the  nonlinear  material 
and  no  intervening  image  forming  components  such  as  lenses  and  mirrors 
exist  between  the  object  and  the  nonlinear  material.  The  upconverter 
aperture  is  limited  by  either  the  crystal  size  or  the  pump  beam  size.  The 
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A.  OPTICAL  SYSTEM  II 


OBJECT  L1  4 


B.  OPTICAL  SYSTEM  I 


FIGURE  10.  IMAGE  UPCONVERTER  SYSTEMS  OF  NON-FOURIER  SPACE 
CONFIGURATION 
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FOV  of  the  system  is  equal  to  the  upconverter  acceptance  angle  0  .  and  the 

lei 


object  size  H  is  limited  to 


H  =  S  tan  0. 


la 


(18) 


From  the  geometry 


tan0imax=  <h  +  a>/2S 


(19) 


Since  object  rays  entering  the  crystal  with  the  incidence  angle  larger  than 
one  half  of  acceptance  angle  will  not  be  effectively  upconverted,  the  upper 
limit  of  0  imax  is  one  half  of  acceptance  angle. 

According  to  the  geometric  optics  analysis,  the  resolvable 

transverse  image  width  of  a  plane  wave  pumped  image  upconverter  operat- 

(2) 

ing  in  the  non-Fourier  space  configuration  is  expressed  by 


Ah  =  D 


1  ~  P 


ni 


tan  0. 

imax 


(20) 


where 


D  =  material  thickness 

/3(=  K./K  =  n.  X  /  n  X.  )  =  Angular  demagnification  factor 

IS  IS  SI 

n.  =  Index  of  refraction  at  m  along  the  phase  match  direction 

0.  =  Angle  between  the  extreme  object  rays  and  the  reference 

imax  J?  J 

axis. 

The  angular  resolution  y  can  be  calculated  using  Equations  (19)  and  (20) 
y  =  tan  (Ah/S)  =  tan  [  D  (1-/3 )  (h  +  a)/2  n^]  (21) 
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Equation,  (21)  indicates  that  resolution  at  the  extremities  of  the  image  is 
poorer  than  that  at  the  center.  The  angular  resolution  at  the  extremities 
and  at  the  center  of  the  image  are  respectively 

fr>  center  *  tan  '*  [Da  (l-fl/2  n.  S2  ]  .  (22) 

i  • 

M  extreme  =ta»'1[D9la(1'-«/2niS],  <22>  ■ 

,  I  t 

It  must  be  remembered  that  equations  (22)  and  (23)  hold  only  for 
the  paraxial  system  (S  >  a).  The  minimum  resolvable  angle  cannot  be 
smaller  than  the  Rayleigh  resolution  angle  of  the  limiting  aperture.  If  the 
crystal  aperture  is  the  limiting  aperture  it  is 

ymin=1-22Va  ■  ' 

The  number  of  resolution  elements  contained  in  the  FOV  is  then 

'  E:Vy 

J 

'  .  I 

The  upper  and  the  lower  bound  of  the  number  of  resolution  elements  for  the 

•  i 

paraxial  system  are: 

aO.  2n,  S2  0 .  2  n.  S  ' 

ia  n  _ _ i  ia  i  (24) 

1.22X.—  Da  (1  -  fij  —  D  (1  -  0) 

=  ,  i 

In  the  optical  system  I,  an  objective  lens  is  used  to  form  an  IR 
image  of  objects  in  the  central  plane  of  nonlinear  crystal  (Figure  10B). 
Under  the  paraxial  assumption  the  extreme  angle  ©imax  is  equal  to  half  the 

angle  subtended  by  the  objective  lens  at  the  nonlinear  crystal 

®imax  =  d/2x£  9ia/2  .  <25>' 

d  =  effective  object  lens  diameter 
x  -  distance  between  the  object  lens  and  the  crystal 


. . . 


Again  0 


is  bounded  by  the  half  acceptance  angle.  The  minimum  resolv 


imax 
able  image  size  is 

Ah'  =  D  (1-/3 )  d/2  n.  <  D  (1-/3)  0.  /2  n. 

'  ^  x  i  '  7  la  i 

The  maximum  number  of  linear  resolution  elements  is 

2ax  m  2am 

R  =  a/Ah'  =  dImR  =  F(MTe  . 

IE 


(26) 


(27) 


For  the  system  imaging  distant  objects  the  distance  x  is  equal  to  the  focal 
length  f i  of  lens  and  the  ratio  x/d  can  be  replaced  by  the  F  number  of 
lens  Lj.  The  upper  bound  of  resolving  power  is  again  limited  by  the  lens 
aperture  diffraction. 


2a  n.  2a  n. 

1.22  X.Pj  >R=  jj  (lJjj)  F1  >  D  . 


(28) 


ia 


Combining  equation  (13)  and  the  right  hand  side  of  equation  (28)  the  lower 

bound  of  the  object  lens  F  number  for  the  distant  object  imaging  system  is, 

D  '1/2 


Fi^ 


X  n 
s  p 


8  X  X.  n.  n 
p  i  l  sm 


(29) 


The  angular  resolution  of  the  optical  system  I  for  upconversion  of  distant 
objects  is 

y=  Ah/fj  =  D  (1-/5)  /2n.  F^  <  D  (1-/3)  0  .Jl  n.  fj  (30) 


The  worst  angular  resolution  is  obtained  when  the  F  number  of  the  object 

lens  is  equal  to  the  lower  bound;  that  is,  F.  =  1/0.  . 

’  1  ia 


The  diffraction  limited  resolution  may  be  obtained  when  the 


F  number  of  the  object  lens  is 


'D(l-g)  * 

.2.44  X.  n. . 
1  1 


1/2 


(31) 
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C.  EFFECTS  OF  CRYSTAL  BIREFRINGENCE 

In  this  paragraph  we  will  consider  effects  of  the  crystal 
birefringence  upon  imaging  properties  of  the  upconverter.  Negligible 
aperture  diffraction  and  plane  wave  pump  are  assumed  in  order  to  isolate 
the  birefringent  effect  from  others.  The  Fourier  optics  formalism  jsed 
in  paragraph  A  is  adopted. 


1.  IMAGE  SPACE  CONFIGURATION 

Omitting  constant  amplitude  and  phase  factors,  the  angular 
spectrum  of  upconverted  field  emanating  out  of  the  crystal  back  surface 


r  /  2 

/  K  K 

6s  <*>  z>  ~  £i  <*’  Sinc  J  [  + 


iffi X  - 

2  x 


1 2  n 

•r  ■  z)\ 


where 


$  (k  ,  z)  = 


z-D/Z  S-D/2  D  D 

K  +  K.  +  2K-.  +  2K 
sa  la  m  f 


The  electric  field  distribution  in  the  upconverted  image  field  becomes 


r 

Es  (Ps  ,  z)  -  d &  e.  (  s,  -s)  sine 


(  K  ,  z)  -  *x  (xg  -  CD/2)  +  Ky  yg 


Equation  (34)  indicates  that  the  optimum  focusing  condition  is 

o 

$  ( x  ,  z)  =  0  for  all  k 

The  solution,  which  defines  the  upconverted  image  plane, 

K.  DIC  r  i  i  i 

?  —  q»  -  __JfL  c  _  sa  1  1  1  1 

Z  S  K  8  2  KL.  +  K  ~  K.  ‘  K 

sa  m  s  1a  sa 


is  identical  to  the  case  of  negligible  birefringence  system  and  defines  the 
image  plane. 

The  line  spread  function  of  a  line  object  parallel  to  the  x  axis 
will  be  considered  first.  The  object  field  distribution  and  its  angular 
spectrum  are 


W  -s)  ~  6(y-yJo) 


£  (  £ ,  -s)  e  y  i0  6  Oc  ) 
1 


The  line  spread  function  is  then 


Ec  (  xo  ,  ye ,  -S')  -  Sine 


Js  v  s  5  Js ’  ~  7  J  UK^K.  y 


k  e 


+i  k  (y  -y.  ) 
yus  ^io ^  d/ 


. . . . .  .  -  . Ml . . . . . . . 


Evaluating  the  integration,  we  obtain 
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+  4 


Mi 


K  D 
P 


-  Sin 


K  K. 
M  i 


v2  r r 


K  D  (ys  yio}  +-?l 
L  P  J 

(37) 


where  Fresnel  functions  S(X)  and  C(y)  are  defined 


S(X)  =  /  sin 


r7T  ,2 
1 


dt 


The  expression  (37)  applies  equally  to  the  system  with  negligible 

birefringence.  This  result  indicates  that  the  image  quality  is  not  degraded 

along  the  y  direction  which  is  normal  to  the  plane  containing  optic  axis  C 

and  the  z  axis  (Figure  11).  The  line  spread  function  is  only  a  function  of 

normal  distance  (y  -  y.  ). 

s  Jio 
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The  line  spread  function  of  a  line  object  parallel  to  the  y  axis  is 
considered  next.  The  object  field  distribution  and  its  angular  spectrum  are 
E.(p.,  -S)  -  6(X.-xio) 

6(  KJ 


ej(*  -S)  e 


X  10  '  y' 


The  line  spread  function  is 
Es(Vxio‘¥-’  ys 


r 

k  2K 

~  &KX  Sine 

D 

2 

x  p  „ 

l2KMKi 

i  k  (x  -x.  - 
xv  s  10  2  ' 


(38) 


Rewriting  above  equation, 


E  (xf -x.  ,  y  ,  -ST)  ^  e 
s  s  io  ’  Js  ’  ' 


.Mx, 

-i  2K  x 


2  I 
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1  K  D 
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where 


(39) 


x 

s 


CD 

2 


Equation  (39)  indicates  that 
image  in  following  ways : 


crystal  birefringence  affects  the  upconverted 


Introducing  change  in  the  field  distribution  from  that  of 
non-birefringent  case 
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•  The  field  distribution  is  a  function  of  (x  -  x.  -  £  D/2) 

s  io 

rather  than  xg  -  x^Q .  This  implies  that  the  x  coordinate 
of  upconverted  image  is  translated  toward  positive  x 
direction  by  an  amount  CD/2,  (Figure  12).  This  linear 
translation  becomes  obvious  if  one  realizes  that  energy 
flow  inside  the  material  is  along  the  anisotropy  direction 
rather  than  along  the  z  direction. 

Line  spread  functions  along  two  orthogonal  directions  expressed  by  equation 
(37)  and  equation  (39)  are  evaluated  for  a  1.06-Mm  pumped  1  cm  long 
proustite  crystal  upconverting  10. 6- urn  objects  into  0. 967-um  images. 

The  results  are  plotted  in  Figure  13.  The  amplitude  of  the  spread  functions 
is  normalized  to  the  peak  amplitude  of  the  line  spread  function  parallel  to 
the  x  axis. 

Resolvable  image  line  width  along  the  y  direction  can  be 
determined  by  Rayleigh  criterion  since  the  sidelobes  of  the  spread  function 
parallel  to  the  x  axis  are  considerably  lower  compared  to  the  mainlobe. 

The  spread  function  parallel  to  the  y  axis  contains  several  nearly 
equal  amplitude  lobes  as  seen  in  Figure  13.  The  amplitude  of  these  lobes  is 
considerably  lower  than  that  of  the  mainlobe  of  the  x  direction  line  spread 
function.  The  width  of  these  lobes  is  also  narrower  than  that  of  the  x  direction 
line  spread  function.  Rayleigh  resolution  criterion  cannot  be  used  to  deter¬ 
mine  resolution  along  the  x  direction,  because  the  adjacent  lobes  are  nearly 
equal  in  amplitude. 
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Resolution  degradation  along  the  x  axis  predicted  from  Figure  12 
has  not  been  observed  in  the  experiments  described  in  paragraph  C  of 
Section  n.  One  of  possible  causes  for  not  observing  resolution  degradation 
is  due  to  the  masking  of  birefringent  effect  by  relatively  larger  diffraction 
line  width  of  the  object  lens. 

More  careful  measurements  of  line  spread  function  parallel  to 
the  y  axis  are  necessary  in  order  to  observe  the  effect  of  crystal 
birefringence. 

2.  FOURIER  SPACE  CONFIGURATION 

It  has  been  shown  in  paragraph  A  that  the  field  distribution  of 
upconverted  image  formed  in  the  lack  focal  plane  of  lens  Lg  (Figure  8)  is 
related  to  the  object  field  distribution  over  the  front  focal  plane  of  lens 
by  repeating  equations  (6),  (7),  and  (8) 


E  (p  }  ^  E.  I 
s  —  s  i 


Ksafl 

Kiaf2 


Sine 
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-l 


AKD 


K  D 
sz 


where  AK  and  K  are  given  by  equations  (7),  and  (8) 
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(6) 


(7) 


(8) 
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From  previous  equations  we  can  conclude : 

*  Magnification  along  two  orthogonal  transverse  directions 
is  not  affected  by  the  crystal  birefringence.  This 
implies  no  distortion. 

•  Different  amplitude  envelope  along  two  orthogonal  trans¬ 
verse  axes.  This  implies  different  acceptance  angle 
for  two  directions.  Acceptance  angle  in  the  yz  plane 

is  given  by  equation  (13).  In  the  xz  plane  it  becomes: 


*pns 

®iax  2  X  n 
s  p 


f 


r2 
£  + 


V 


2npXs  Xj 

X  n  n. 
p  x  1 


1/2 


-e 


•  No  degradation  in  resolving  power,  since  the  factor 

AKD 


Sine 


image. 


in  equation  (6)  only  affects  the  envelope  of 


(40) 


D.  POINT  SPREAD  FUNCTION  OF  A  GAUSSIAN  BEAM  PUMPED 
FOURIER  SPACE  CONFIGURATION  SYSTEM _ 

The  effect  of  pump  beam  amplitude  distribution  on  the 
imaging  characteristics  of  the  Fourier  space  configuration  image  upconverter 
is  considered.  Diffraction  due  to  the  lens  apertures  and  the  crystal  birefrin¬ 
gence  are  neglected.  It  is  also  assumed  that  the  physical  size  of  the  non¬ 
linear  material  is  large  compared  to  the  pump  beam  size.  The  optical 
arrangement  is  shown  in  Figure  8. 


50 


I  -i 


1  I 


M 


The  pump  beam  is  assumed  to  be  in  the  TEM  mode.  The 

oo 

Gaussian  window  at  the  entrance  surface  of  the  nonlinear  crystal  is 


expressed  by* 


2/  2 

E  (p,  -D/2)  ^  e  /a 


The  angular  spectrum  of  pump  field  distribution  is 


ep(  Ep,  -D/2)  -  e 


-■  *,■/* 


The  point  object  in  the  front  focal  plane  of  the  object  lens  is  expressed  by 


E.(o.,  -2frD/2)  - 


The  angular  spectrum  in  the  IR  object  fields  in  the  entrance  plane  of 
the  crystal  is,  using  equation  (4), 

£.(  £., -D/2)  -  E.(-L  £./K.  ,  -2L-D/2)  <~5(£.  +K  p.  /L)  (45) 

i'  i’  l  1  l  1a  1  i  ia-10  1 


The  sum- frequency  polarization  P  is  obtained  by  taking  a  convolution  of 

b 

£.  (  k.  ,  -D/2)  and  E  {  k  ,  -D/2). 
i  i*  P  P 


♦Constant  amplitude  and  phase  factors  are  again  omitted  in  the  following 
analysis. 
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V  Ss,-D/2 )~J  e.(  if.,- -D/2).  £p(  *s-*p,  -D/2)  d  *. 

-a2(  /c  +K.  p.  /f,)2/4 
v  ~s  ia-io'  i  / 


The  angular  spectrum  of  upconverted  image  fields  at  the?  exit  surface  of  the 


crystal  is 


-a2(  S. 

es (  “s ’  +D/2)  ~  e 


+K.  p../L)2/4  ^AimN  -1d(ko 

m-ld  1  '  Sinc^PJe  ^sz  2  ' 


where  £K  and  K  are  mismatch  factors  and  the  z  component  of  wave 

i  SZ  !  S  ; 

vector  in  the  material  for  the  angular  spectrum  with  the  angular  . 
frequency  k  g .  , 


Field  distribution  of  the  upconverted  image  in  the  back  focal  plane  of  the 


imaging  lens  is,  by  equation  (5) 


Es  <_os  ,  2  f2  +D/2)  ~  6s  (K6a^6/f2 ,  +D/2) 


0  f  K  K.  ,  ' 
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The  symbol 


*  =“C-°s 


indicates  again  that  AK  and  K  are 

sz 


s  *2 

evaluated  by  replacing  £0  with  K.  /f„  as  shown  in  equations  (7)  and  (8). 

s  sa  —  s  z 

Equation  (48)  is  the  point  spread  function  of  Gaussian  beam 

pumped  Fourier  space  configuration  image  upconverter.  The  upconverted 

image  is  no  longer  a  point.  The  sine  function  again  limits  the  FOV  of  the 

system  and  describes  the  amplitude  envelope  of 1  he  image  plane  in  the 

similar  way  as  the  plane  wave  pumped  case.  The  Gaussian  amplitude  factor 

o 

,2 


aK. 


la 
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1  sa 


2  L 


0_  +  o.l 


f2Kia  ~s  '~ioy 


indicates  that  transverse  magnification  factor  M,  and  the  image  spot  size  are 
given  by 


Mt  =-s/-i  =  “f2K 


■J- 


Ksa 


(49) 


Ws  =2f2/aKsa 


(50) 


Equation  (49)  is  identical  to  equation  (11).  It  implies  that  the  Gaussian 
amplitude  distribution  does  not  introduce  distortion.  Combining  equations 
(49)  and  (50)  and  dividing  the  object  spot  size  by  the  focal  length  of  the  object 
lens  we  obtain  expressions  for  the  resolvable  IR  object  spot  size  and  angle; 


W.  *  X .  f.  /  7ra 

l  i  1; 


(51) 
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Since  Rayleigh  resolution  criteria  cannot  be  applied  to  the  present  case 
the  1/e  amplitude  is  arbitrarily  chosen  as  the  resolution  limit. 


From  above  discussion  one  may  conclude: 

•  FOV  and  the  intensity  envelope  of  image  plane  may  not 
be  affected  by  the  pump  beam  intensity  distribution. 

•  The  optical  transfer  function  of  the  Fourier  space 
configuration  system  depends  heavily  on  the  pump  beam 
intensity  distribution.  A  multimode  pump  beam  (diverging 
pump  beam),  therefore,  degrades  resolution  of  the  Fourier 
space  configuration  upconverter. 


IV.  UPCONVERSION  EFFICIENCY 


Conversion  efficiency  of  the  upconversion  process  has  been 
measured  and  compared  to  a  theoretical  calculation  based  on  the  plane  wave 
theory  by  several  authors^5  Measured  efficiency,  the  conversion 

efficiency  formula  used  for  checking  their  results,  ?nd  parameter  values 
such  as  SHG  d  coefficients  are  inconsistent  among  -nvestigators.  In  this 
section,  an  attempt  will  be  made  to  clarify  some  of  the  confusion  by  con¬ 
sidering: 

•  Proustite  SHG  d  coefficients 

•  Experimental  measurement  of  upconversion  efficiency 
under  controlled  experimental  conditions 

A.  PROUSTITE  SHG  d  COEFFICIENTS 

The  only  reported  measurement  of  SHG  d,^  coefficients  of 

proustite  was  performed  by  K.  H.  Hulme  et  al  by  comparing  SHG  of  1. 152  Mm 

(12) 

with  proustite  to  SHG  with  LiNbO^  along  the  crystallographic  X-axis'  ' .  The 

ratio  ldggl  /  Id^^l  of  proustite  was  determined  by  Hulme  et  al  by  the  fre- 

(12) 

quency  summing  technique  .  Their  measurement  did  not  determine  the 
absolute  signs  of  the  coefficients.  Their  measured  results  are 

(53) 
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3 


,P/dQ1P  =1.6  ±0.1 

5  /  ul 


The  inconsistency  in  proustite  dg^  is  resulted  from  the  spread 
in  LiNbOg  dg^  coefficient  reported  in  the  literature  as  shown  below. 


G.  Boyd 


(13) 


d31,LN)  =  (10.6±1.0)d36(KDP) 


R.  Miller  and  Savage 


(14) 


d31(L!''  =  (11.9  ±1.7)  d36(KDP) 


Bjorkhc'm*15’ 16) 


d31(LN)  =  (10.9  ±1.7)  d3f.(KDP> 


(17) 

Byer  and  Harris  ' 


d31<LN)  =  11.0  10.6  d36(KDP> 


Careful  review  of  the  above  reported  measurements  indicates  that  Bjorkholm’ 
measurement  seems  more  accurate  than  others. 

Using  Jerphagnon  and  Kurtz anc}  Francois^^^  measurement 
on  dgg  of  KDP  crystal, 


A  (KDP) 
a36 


=  1.67  x  l(f9  esu 
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and  Bjorkholm’s  d, 


one  finds 


(LN) 


31 


id*. 

i  3i 


(P) 


-  28.4  d36(KD)  =  42.6  x  10“9  esu 


(54) 


(P) 


22 


=  45.  5  d36(KDP)  =  68. 2  x  10“9 


esu 


(12) 

Hulme  '  et  al  gave  estimated  d  coefficients  from  their  measurement  as 


*22 


(P) 


*31 


(P) 


=  50  d36(KDP)  =  75  x  10"9  esu 


=  30  d36(KDP)  =  45  x  10“9  esu 


(55) 


The  reason  for  this  modification  from  experimental  results  was  not  given  by 
Hulme. 

Boggett^9^  et  al  reported  the  effective  d  coefficient  of  Type  I 

phase  matched  proustite  by  measuring  doubled  10.6  4m  output  power.  Their 

results  for  0  =  22. 5°. 

m 


60  d36 


(KDP) 


is  higher  than  the  calculated  value  based  on  Hulme's  estimated  d  coefficients 


(KDP) 

=  30  d36v  ' .  If  one  solves  the  equation 


deff  =  d22  008  9m  +  d31  sLn  8m 
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which  applies  to  Type  I  phase  matching,  assuming  ^32^31  =  1*  6  for 

(P)  (imp) 

proustite,  one  finds  d^'  ;  =32.2  d^g'  '  from  the  Boggett  measurement. 
Boggett’s  result  is  very  close  to  the  Hulme’s  estimated  d^^-  value. 

_Q  _Q 

The  values  d ^  =  68. 2  x  10  esu  and  d^^  =  42. 6  x  10  esu  will 
be  used  for  following  system  performance  calculation,  although  Hulme’s 
estimated  d  coefficients  are  widely  adapted  by  other  workers. 

The  choice  is  dictated  by  the  conversion  efficiency  measurements 
performed  at  AIL.  The  result  of  the  most  recent  conversion  efficiency  mea¬ 
surement  is  presented  in  paragraph  B. 


B.  CONVERSION  EFFICIENCY  MEASUREMENTS 

The  main  purpose  of  the  experiment  is  to  assess  the  achievable 
conversion  efficiency  of  a  practical  image  upconverter  as  compared  to  the 
theoretical  estimate.  The  experimental  conditions  for  efficiency  measure¬ 
ments  must  be  carefully  controlled  in  order  that  the  experimental  results 
can  be  meaningfully  compared  to  the  theory. 

Conversion  efficiency  of  the  experimental  model  was  measured 
with  TEM  mode  signal  beam  and  TEM  mode  pump  beam.  The  uocon- 
verted  signal  was  detected  by  a  low  noise  SI  photomultiplier  with  accurately 
calibrated  quantum  efficiency.  Optical  efficiency  of  optical  components  such 
as  interference  filter,  optical  windows,  and  the  beam  combining  dichroic 
mirror  were  carefully  calibrated.  The  pump  beam  was  chopped  with  a  duty 
factor  of  one  quarter  to  prevent  excessive  crystal  surface  oxidation. 
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The  input/output  linearity  was  checked  so  that  an  efficiency  mea¬ 
surement  was  taken  within  the  dynamic  range  of  the  system.  Figure  14  shows 
the  plot  of  input  C09  power  versus  the  upconverted  signal  induced  pho4:^ elec¬ 
tron  at  the  photocathode.  The  observed  power  upconversion  efficiency  is 
calculated  by 


N  hC 

_  s _ 

\  =  T)  d  P.  X 
'p  lr  s 

where 

Ng  =  photoelectron  flux 

17  =  quantum  efficiency  of  the  photomultiplier 

r 

h  =  Planck's  Constant 

C  =  velocity  of  light 

X  =  wavelength  of  upconverted  signal 
s 

P.^  =  IR  input  power 
d  =  duty  factor  of  light  chopper 
Using  the  following  values: 

d  =  0. 25 

*Y~=  2.05  x  10  Joules  for  X0  =  0.967  nm 
s  s 

r?p  =  9  x  10“4 

Ns  =  4x  105  Hz 

P.  =0.09  W 
it 


(56) 
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-9 

the  measured  system  peak  power  conversion  efficiency  is  4.05  n  10  which 

is  approximately  one  third  of  the  theoretical  estimated  efficiency  of 
-*8 

1.45  x  10  calculated  by  the  expression 


104.44  d  J  l2  P 
eff 


2 

1  - - =5 - ^  -  2. 162  X  lO^-fe'f 

"i  ns  np  Xs  <Ap  +  Ai>  '  p  V 


(57) 


where 

P  =  pump  power  in  watts 

Jr 

l  =  interaction  length  in  cm 
n^  =  indices  of  refraction  at  (l  =  i,  p,  s) 

Ap  =  pump  beam  area  in  cm^ 

O 

A.  =  IR  beam  area  in  cm 

t  =  optical  transmission  of  the  experimental  model 

Parameter  values  used  in  the  efficiency  estimate  are: 

-9  -9 

d  ,,  =  d„„  cos  0  =68  (cos  20°)  x  10  esu  =  64  x  10  esu 

eff  ll  m 

n.  =  2.68 

np  =  2. 817 

n  =2.804 
s 

l-  0.6  cm 
X  =  0.967  4m 

b 

o 

A  =  1.13  cm  (beam  diameter  =  1.2  cm) 

hr 
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► 


2 

=  0.783  cm  (beam  diameter  =  1  cm) 
=  4.8W 
=  7.44  x  10'3 


The  most  probable  cause  for  the  discrepancy  between  the  mea¬ 
sured  efficiency  and  the  prediction  based  on  equation  (57)  is  the  phase  mis¬ 
match  due  to  the  beam  divergence.  Equation  (57)  was  derived  for  the  case  of 
perfectly  phase  matched  interaction  of  two  near -field  Gaussian  beams.  Other 
factors  affecting  efficiency  reduction  are  the  beam  walk-off,  the  crystal 
birefringence,  and  the  crystal  absorption. 

System  conversion  efficiency  under  the  experimental  condition 
that  nearly  realizes  the  analytical  model  was  measured.  The  measured  con¬ 
version  efficiency  is  approximately  a  factor  of  1.75  less  than  the  theoretical 
prediction.  In  this  measurement  the  10.6  jim  beam  was  spatially  filtered  by 
a  small  pinhole  placed  at  beam  center.  The  filtered  10.6  (im  wave  front  is 
very  close  to  plane  wavefront.  Beam  divergence  effect  is  minimized  by 
aligning  the  optics  so  that  the  interaction  of  DR  waves  with  the  pump  waves 
is  limited  to  the  center  of  the  pump  beam. 

Generally,  the  conversion  efficiency  formula  used  for  the  system 

calculation  is  different  from  equation  (57).  The  generally  accepted  formula 

.s(21, 22) 


h  = 


52.22  d  3  l2  t 

eff  op  p 

o~ 

n.  n  n  X  “  A 
i  s  p  s  p 


(58) 
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which  is  obtained  under  the  assumption  of  uniform  intensity  pump  and 
A  >>  A.  «  A  . 

p  is 

The  measured  conversion  efficiency  with  a  small  IR  signal  beam 
described  here  is  1.75  times  lower  than  the  theoretical  estimate  based  on 
equation  (57),  but  is  within  15  percent  of  the  estimate  based  on  equation  (58). 
It  must,  however,  be  noted  that  the  measured  conversion  efficiency  is  the 
peak  efficiency,  since  the  interaction  takes  place  at  the  center  of  pump  beam, 
where  the  power  density  is  the  highest.  Reduction  in  the  efficiency  at  the 
skirt  of  the  pump  beam  is  expected. 

From  the  experiments  described  and  considerations  made,  it 
seems  reasonable  to  conclude  that  power  conversion  efficiency  estimated  by 
equation  (58)  can  be  practically  achieved  at  the  center  of  FOV  of  the 
upconverter. 


V.  PLANS  FOR  EXPERIMENTS  AND  ANALYTICAL  STUDIES  FOR  THE 

REMAINDER  OF  THE  PROGRAM 


The  experimental  effort  for  the  remainder  of  the  program  will  be 


directed  toward: 


Optimization  of  a  pulse  gated  image  upconverter  system 
and  measurement  of  its  operating  characteristics  such 
as  conversion  efficiency,  resolving  power,  etc. 

Crystal  damage  experiments  due  to  a  high  power  pulse 
pump  with  moderate  repetition  rate. 

Measurement  of  thermal  conductivity  proustite. 

Image  upconversion  experiments  with  diffuse  objects. 

Experimental  verification  of  the  crystal  birefringent 
effect  on  the  line  spread  function 

The  following  areas  will  be  studied  analytically: 

Effect  of  pump  beam  intensity  distribution  on  the  image 
intensity  envelope  for  the  non- Fourier  space  configured 
upconverter. 

Accurate  conversion  efficiency  calculation  of  the  image 
space  configuration  image  upconverter. 

System  design  calculation  of  the  range  gated  image  up¬ 
converter. 

Estimation  of  thermal  detuning  due  to  pump  absorption. 
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